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Abstract 

Using the Behe detector operating at the KEKB e'^e~ storage ring, we have measured the mean 
muhiphcity and the momentum spectrum of neutral pions from the decays of the T(4S') resonance. 
We measure a mean of 4.70 ± 0.04 ± 0.22 neutral pions per T(4S') decay. 
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This analysis presents the first results on the average multiplicity of neutral pions and 
their momentum spectrum in B meson decays at the T(4S') resonance using the Belle de- 
tector [0] at KEKB 0. The measurement is based on the sample of 3.4 x 10^ B meson pairs 
collected by the Belle detector at a center-of-mass (CM) energy of y/s = 10.58 GeV. 

The particle composition of hadronic final states in e^e~ annihilation and the measure- 
ment of inclusive particle production rates have been important subjects for various energy 
regions. Previous studies of particle composition at the T(45') resonance include measure- 
ments of charged particles (tt^ Q and K"^ @5 0), V mesons [0], and vector mesons 
For typical B meson decays, the bulk of the neutral energy is carried by neutral pions. Mea- 
surements of inclusive spectra contain information on B meson decay mechanisms, especially 
in the high momentum region where important rare decays may become detectable. If the 
details of these inclusive spectra are known more precisely, they will allow better estimation 
of backgrounds and better modeling of B meson decay. 

The Belle detector (see Fig. ||), described in detail in Ref. consists of a silicon vertex 
detector (SVD) 0, central drift chamber (CDC) [|lTJ, aerogel Cerenkov counter (ACC) [ll2|, 
time of flight /trigger scintillation counter (TOF/TSC) |jl3[, Csl electromagnetic calorimeter 
(ECL) fl^ and K^/muon detector (KLM) [l^. The SVD measures the precise position of 
decay vertices. It consists of three layers of double-sided silicon strip detectors (DSSD) in 
a barrel-only design and covers 86% of the solid angle. The layer radii are 3.0, 4.5, and 6.0 
cm. Charged tracks are reconstructed primarily by the CDC that covers the 17° < 6 <150° 
polar angular region. It consists of 50 cylindrical layers of drift cells organized into 11 
super- layers (axial or small-angle-stereo), each containing between three and six layers. A 
low Z gas mixture (50% He, 50% C2H6) is used to minimize multiple- Coulomb scattering. 
The inner and outer radii of the CDC are 9 cm and 86 cm, respectively. The solenoidal 
magnetic field of 1.5 Tesla is chosen to optimize momentum resolution without sacrificing 
reconstruction efficiency for low momentum tracks. Kaon identification (KID) is provided by 
specific ionization {dE/dx) measurements in the CDC, Cerenkov threshold measurements in 
the ACC, and the cylindrical TOF scintillator barrel. The ECL is made of finely segmented 
CsI(T£) crystals 30 cm in length. The cross section of one counter is approximately 55 x 55 
mm^ at the front surface. The ECL crystals cover the 12° < 6 < 157° angular region. The 
inner radius of the barrel is 1.25 m, while the annular endcaps are placed at +2.0 m and 
— 1.0 m along the beam line from the interaction point. The calibration of the calorimeter 
is performed using cosmic rays and Bhabha events. The KLM consists of alternating layers 
of glass resistive plate counters and 4.7 cm thick iron plates. 

The data samples used in this analysis correspond to 3.2 fb~^ of integrated luminosity 
taken at the T(4S') resonance and 0.6 fb~^ taken at a CM energy 60 MeV below the reso- 
nance; the latter was used to subtract underlying continuum background. The integrated 
luminosity was determined from the number of Bhabha events for which we require both 
electron and positron in the region of 46.7° < 6* < 145.7° in the center-of-mass frame. 

Hadronic events are selected based on charged track information from the CDC and clus- 
ter information from the ECL. We require at least three charged tracks, that the energy 
sum in the calorimeter be between 10% and 80% of y/s, and that the charged track mo- 
mentum be balanced in the z direction. This removes the majority of two photon, radiative 
Bhabha, and r+r~ events where both r's decay to leptons. Radiative Bhabha events with 
one electron outside of the ECL acceptance are removed by requiring at least one large-angle 
cluster in the ECL and requiring that the average cluster energy be below 1 GeV. Higher 
multiplicity t~^t~ events are removed if the charged and neutral energy sums in the event are 
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consistent with a r pair event and if the reconstructed invariant mass of the particles found 
in each of the two hemispheres perpendicular to the event thrust axis falls below the r mass. 
Beam- gas and beam- wall backgrounds are removed by reconstructing the primary vertex of 
the event and requiring it to be consistent with the known location of the interaction point. 
The selection is 99% efficient for BB events and approximately 87% efficient for continuum 
events. The efficiency for the hadronic event selection was found by a GEANT-based |T^ 



Monte Carlo simulation program. To suppress continuum, we require that the ratio R2 of 
second to zeroth Fox- Wolfram moments ||T^, determined using charged tracks and neutral 
clusters, be less than 0.5. 

Photons are reconstructed from neutral clusters in the ECL that have a lateral shape 
consistent with that of an electromagnetic shower. The energy resolution was measured 
to be aE/E = 0.066%/E © © 1.34% {E in GeV) from beam tests H. To 

keep the combinatorial background at a reasonable level, only photons in the central barrel 
region (35° < 9 <120°) with E^ > 30MeV are used in this analysis; the endcap regions 
have worse energy resolution due to more intervening material and higher beam-associated 
background. 

For each 100 MeV/c momentum bin in the CM momentum range to 3 GeV/c, the 77 
invariant mass distribution is fit to an asymmetric (symmetric above 2 GeV/c) Gaussian 
(i.e., a Gaussian with different widths on either side of the mean) for the signal plus a 
polynomial for the combinatorial background to extract the 7r° yields. Fig. ^ shows typical 
mass spectra obtained from the on-resonance data. For the asymmetric Gaussians, the 
mass resolution is defined as the mean of the left- and right-hand sigmas. An average mass 
resolution of 5 MeV/c^ is obtained, dominated by energy resolution at low momenta or by 
angular resolution at high momenta. The mass peak is shifted slightly from the established 



value |T8| because of the asymmetric energy response of the calorimeter due to shower 
leakage. The observed mass peak position and resolution are consistent with Monte Carlo 
expectations, as shown in Fig. |^. 

To extract the 7r° momentum spectrum from T(4S') decays, the underlying continuum 
in the on-resonance data is subtracted bin-by-bin using off-resonance data. The inclusive 
spectrum is calculated using 



where Nh is the number of produced BB events, yon and Yos are the background-subtracted 
7r° yields obtained from on- and off-resonance data fits, a = {Con/ i^os) • (sofr/son) is the 
on-off scaling factor, and e is the product of acceptance and detection efficiency for each 
momentum bin. The average vr*^ multiplicity is obtained by summing the data from the 
measured individual momentum bins as shown in Table |. 

The TT*^ acceptance and detection efficiency are determined from Monte Carlo simulations 
of BB decays (equal proportions of charged and neutral B mesons) and continuum processes. 
The 77 invariant mass distributions are fit with the same functions as used in the real data 
analysis. The product of acceptance and detection efficiency are defined as the ratio of the 
fitted TT^ yield to the generated count. Efficiencies for the high momentum 7r°'s above the 
kinematic limit for BB decays are deduced from continuum Monte Carlo normalized in the 
2.2-2.3 GeV/c bin, and scaled efficiencies are used that account for the different acceptances. 

Possible sources of systematic uncertainties and their effect on the inclusive vr" mean 
multiplicity measurement are summarized in Table These are discussed next. 
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The uncertainties in the number of BB events and the hadronic event selection efficiency 
are estimated to be 1% and 1.1%, respectively, or 1.5% combined. 

The effect of uncertainty in the relative luminosity ratio, used to subtract the continuum 
background from the on-resonance data, was studied by varying the size of the continuum 
subtraction by 1%, and we find a 1.5% uncertainty in the mean multiplicity. 

Due to the CM energy difference between on- and off-resonance data, the characteristics 
of continuum events may not match, leading to a bias in the continuum subtraction from the 
on-resonance data. In particular, the event shape at lower CM energy is slightly less jet-like 
so that more continuum events survive the common R2 cut (resulting in an oversubtraction 
from the on-resonance data), and the particle momenta scale with CM energy. Hadronic 
event selection and R2 cut bias from the on-off energy difference were determined by using 
a continuum Monte Carlo [|l^] sample generated at the CM energy of the off-resonance 
data; we estimate a 0.4% uncertainty in the multiplicity due to this effect. Particle 
momentum scaling with CM energy was studied by comparing the multiplicity with and 
without momentum scaling; this leads to an uncertainty of 0.3%. 

The uncertainties in the tt*^ detection efficiency due to the application of the shower 
transverse shape cut, the charged-track veto, photon energy smearing, and the ECL's mod- 
est non-linear energy response correction, are estimated to be 0.4%, 1.6%, 2%, and 0.6%, 
respectively (determined by comparing the yields with and without each of these criteria). 

The decay angular distribution for the photons in the 7r° rest frame is expected to be 
isotropic, but the detected distribution of the decay angle 9d — defined as the angle between 
the photon momentum in the 7r° rest frame and the tt^ momentum in the laboratory frame — 
shows an energy dependence due to the 7 energy cut. Inaccuracy in the Monte Carlo 
simulation could manifest itself as an anisotropy in the efficiency-corrected decay angle 
distribution. A comparison between Monte Carlo and data for | cos6'd| > 0.5 and | cos 6*^1 < 
0.5 shows no discrepancy and is used to estimate a systematic uncertainty of 0.3%. 

The 77 invariant mass distribution for each CM momentum bin is fit to the following 
functions in the mass window 80-180 MeV/c^: 

• p*^ < 1.0 GeV/c: 4th order polynomial plus asymmetric Gaussian; 

• 1.0 < p*^ < 2.0 GeV/c: 2nd order polynomial plus asymmetric Gaussian; 

• 2.0 < p*^ < 3.0 GeV/c: 2nd order polynomial plus Gaussian. 

The fitting range and the order of the polynomial are chosen to minimize statistical fluctu- 
ation. By changing the order of the polynomial, the relative yield variation between data 
and Monte Carlo is taken as the systematic uncertainty in the background modeling; a 2.8% 
uncertainty was deduced. The weighted average of the fit errors in the efficiency estimation, 
(5e)^ = J2i ^li'Y^ I Yl,i ! leads to 1% uncertainty in the 7r° multiplicity. Combining these 
two numbers gives a 3% systematic uncertainty in the 7r° mean multiplicity measurement 
due to the fit procedure. 

The measured 7r° momentum spectrum is compared to the BB Monte Carlo |20| predic- 



tion in Fig. ^ Systematic differences in the intermediate momentum region may be caused 
by overestimation of 6 c processes in the Monte Carlo event generation: 100% 6 ^ c is 
assumed for generic BB decays in our Monte Carlo. In the high momentum region above 
the b ^ c kinematic end point, we searched for a 7r° excess from charmless B decays. We 
find Nexcess = 410 ±724 events in the momentum interval 2.4-2.7 GeV/c, which corresponds 
to a partial branching fraction B{B 7r°X; p > 2.4GeV/c) = (2 ±3.5 ±6.7) x 10""^, or less 
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than 1.4 X 10~^ at 90% confidence level. However, the statistical precision here is limited 
and the fluctuations are dominated by the off-resonance data used in the subtraction. 

In conclusion, using 3.2 fb~^ of data accumulated at the T(45') resonance by the Belle 
detector, we have measured the inclusive spectrum of neutral pions from T(45') decays. By 
summing the measured momentum bins, the mean multiplicity of neutral pions from T(4S') 
decays has been determined to be (n^o) = 4.70 ± 0.04 ± 0.22, corresponding to an inclusive 
branching fractional] of B{B — > tt^X) = (235 ± 2 ± 11)%, where the first error is statistical 
and the second is systematic. 
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FIG. 1: Isometric cutaway view of the Belle detector. 
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FIG. 2: 77 invariant mass distributions for the T(4S') resonance data, (a) p*^ = 0.0-3.0 GeV/c; (b) 
0.0-1.0 GeV/c; (c) 1.0-2.0 GeV/c; (d) 2.0-3.0 GeV/c. An average mass resolution of 5 MeV/c^ was 
obtained. The smooth curve in each plot is a fit to the data using an asymmetric (or symmetric) 
Gaussian plus a polynomial. 
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FIG. 3: Mass peak position and resolution obtained from ofF-resonance data and Monte Carlo 
expectations. (a),(c) data; (b),(d) continuum Monte Carlo. The fluctuations near 2 GeV/c are 
caused by the change in the signal fitting function from an asymmetric to a symmetric Gaussian. 
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FIG. 4: Measured inclusive ir^ momentum spectrum in T(45') decays, with the high momentum 
range shown in the inset. The histograms show the Monte Carlo prediction. The error bars are 
statistical, while the shaded band indicates one standard deviation systematic uncertainties. 
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TABLES 
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TABLE I: Measured inclusive 7r° spectrum from T{AS) decays, using 3.2 fb on-resonance and 
0.6 fb~^ off-resonance data. 
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1 

it 


U.UZD 


1 


.8-1.9 


10759. 


Q 
.O 


lb 


805. 


Q 

,o 


37.0 


0.084 


1 

it 


U.UUo 


1 

it 


U.UZZ 


1 


.9-2.0 


7842, 


7 
. ( 


lb 


746. 


Q 

,o 


36.7 


0.061 


1 

it 


U.UUo 


1 

it 


U.Uio 


2, 


.0-2.1 


7241, 


.6 


lb 


667, 


Q 

,o 


35.2 


0.059 


1 

it 


U.UUo 


1 

it 


n ni 7 
U.Ui ( 


2, 


.1-2.2 


4945, 


.7 


lb 


586. 


,1 


32.4 


0.044 


lb 


0.005 


lb 


0.014 


2, 


.2-2.3 


2956, 


.1 


lb 


533. 


,1 


31.2 


0.027 


lb 


0.005 


lb 


0.011 


2, 


.3-2.4 


-345, 


.5 


lb 


505. 


,8 


30.7 


-0.003 


lb 


0.005 


lb 


0.011 


2, 


.4-2.5 


163, 


.3 


lb 


443. 


,3 


29.7 


0.002 


lb 


0.004 


lb 


0.009 


2, 


.5-2.6 


488, 


.3 


lb 


415. 


,3 


28.4 


0.005 


lb 


0.004 


lb 


0.007 


2, 


.6-2.7 


-242, 


.0 


lb 


394. 


,3 


27.2 


-0.003 


lb 


0.004 


lb 


0.007 


2, 


.7-2.8 


189, 


.8 


lb 


347. 


,3 


25.5 


0.002 


lb 


0.004 


lb 


0.006 


2, 


.8-2.9 


-206, 


.9 


lb 


314. 


,9 


23.7 


-0.003 


lb 


0.004 


lb 


0.005 


2, 


.9-3.0 


-59, 


.5 


lb 


278. 


,8 


21.6 


-0.001 


lb 


0.004 


lb 


0.005 
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TABLE II: Sources of systematic uncertainty in the inclusive tt multiplicity measurement. 



Source Effect on (n^o) (%) 



Overall normalization 


1.5 


= ^on - (f^) (If) • i^off 


1.5 




0.4 


OFF E* ■ .[^ 


0.3 


non-linearity 


0.6 


E-^ smearing 


2.0 


Shower shape 


0.4 


Track match 


1.6 


Fit procedure 


3.0 


Energy dependence of efficiency 


0.3 


Total (added in quadrature) 


4.6 
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